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1. INTRODUCTION 

The practitioner has to think about (at least) two obstacles, when he wants to 
solve a Hamiltonian field theory problem in front form dynamics. One is the 
problem that the theory is not manifestly rotational invariant. The other is the 
question, how one should construct effective theories at all in this framework. 
Of course, the first problem exists analogously in instant form dynamics, where 
the theory is not manifestly Lorentz invariant. Connected to this problem is 
the fact that any symmetry which is not manifest will break down immediately 
with the slightest approximation in a Hamiltonian field theory. There is broad 
agreement that the construction of effective theories is inevitable, if one wants 
to describe the low energy region, i.e. the bound states, of a strongly coupled 
theory. A variety of formalisms were suggested concerning this topic [H |^ . 

In this article I shall present quantitative investigations concerning the above 
two problems with a specified model, namely positronium as a QED(3+1) 
bound state. The model is represented by the energy diagrams of Fig. |^. To 
be more specific, an effective theory was constructed following the work of 
Pauli 1^. The effective matrix elements were calculated analytically, and the 

(*) Invited talk at the International Workshop “New Non-Perturbative Methods 
and Quantization on the Light Cone”, Les Houches (France), February 24-March 7, 
1997; to appear in the proceedings. 
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Fig. 1. — The graphs of the positronium model. Effective photon lines are labeled 
by hashed rectangles. The graphs of the annihilation interaction are at the bottom 
line. 


emerging integral equation was put on the computer to numerically solve for 
the eigenvalues and eigenfunctions of the Hamiltonian. 

Proceeding in this manner, one can answer the following questions quantita¬ 
tively: Are the results of front form theory rotationally invariant, i.e. do states 
of different Jz form multiplets? How good is the underlying effective theory 
concerning its cutoff dependence? Is the multiplet structure of the spectrum 
reproduced correctly? And: what is the numerical evidence for the degeneracy 
of corresponding states? 


2. METHOD 


A complete description of the applied method is surely beyond the scope of 
this article. I review shortly its major steps. 

To get a finite dimensional Hamiltonian out of the infinite dimensional canon¬ 
ical Hamiltonian, we introduce a cutoff A on the kinetic energies. We then 
(formally) map an x onto a 2 x 2 block matrix by subsequent projections 
(Bloch-Feshbach formalism) and call these sectors P-space and Q-space. Note 
that the Q-space is already an effective sector by construction: 

TT _ _L 

’ 


where 


is a piece of the original Hamiltonian and ATg® are the interactions 
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generated by the projections. The Hamiltonian reads 

The final projection maps the Q-space onto the P-space. One solves for the 
state Q\'il}) = leey) with a resolvent involving the redundant parameter uj 


G{uj) := {Q\aj — PlcIQ) 


-1 


and obtains the nonlinear equation 

with the effective Hamiltonian (P|V') = |ee)) 

Hl%iu) := PPlc-P + PPlcQ(^ - Hi^c)~^QHi^cP- 


The fixing of the redundant parameter works as follows. We divide the sector 
Hamiltonian in Q-space into Hq = Mq+Vq, and the interaction into a diagonal 
part (Vq) and a non-diagonal part 6Vq 

= (^q) + ^^Q- 


With the definition T* := oj—{Vq) we expand the resolvent around the diagonal 
interaction (Vq) 


1 


1 


1 




1 


CO -Hq T* - T* - ^T* -M^- SVq 

1 


T* - M^' 


( 1 ) 


The approximation in the last step is to consider the first term only. It has 
severe consequences: a collinear singularity occurs, proportional to 


1 A(x,k±,x' ,k'j;T*) 


with V{x,x';T*)— {x - x'\{T* - Ml^). 


V \x — x'\ 

We can calculate A{x, kx,x', T*) with P-space graphs 


where := (fc' — and ?g := {kg — k'g)^ are the momentum transfers. 
One determines the parameter uj by demanding that the collinear singularity, 
induced by the approximation, Eq. (|^), vanishes 

A{x,kj_,x',k'j_;T*) = 0, '\/x,x ,k±,k'^. 
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This allows for the calculation of an explicit expression for T* 


T*{x, k±;x',k'j_) 


1 / mj + k± 

2 I x{l — x) 


rri^f + k'^\ 
x'{l — x')j ' 


The interpretation is that T* is an approximation of the summed interactions 
of the higher Fock states. 

Finally, the procedure yields an integral equation. The effective Hamiltonian 
operates only in P-space, and the continuum version of the eigenvalue problem 
is 

/ mi + \ 

0 = f j ^/>„(a;,fc_L;Ai,A2) 

f dx'cPk'^ j^(4,Ae)j^(/g, Ag) 

IBtt^ ^xx'{l — a;)(l - x') lele.ii 

Ai ,A2 

xtpn{x', k'j_; A^, A 2 ), 


with an 

states. 

Lorentz 


integration domain D defined by a cutoff on the kinetic energy of the 
One observes that the effective interaction is gauge invariant and a 
scalar 


UeB ■ = 


Ae)j/i(^e) Ae) 


/M; 


3 . RESULTS 

As pointed out in the last section, the goal is to construct an effective Hamil¬ 
tonian describing the positronium bound states. We showed how to find a 
reasonable ansatz for contributions of higher Fock states. To find out about 
the rotational symmetry of the theory, one has to observe that the rotation 
operator J 3 around the z axis is kinematic. Consequently, states can be clas¬ 
sified according to the corresponding quantum number J^. We have shown 
in Refs. 0[|| that closed expressions for the matrix elements of the effective 
Hamiltonian can be calculated, even for arbitrary Jz- This enables one to solve 
for the positronium eigenvalues in each sector of Jz separately and to compare 
the results to find out about possible degeneracies. A computer code for solv¬ 
ing these eigenvalue problems, i.e. integral equations, has been generated with 
correct Coulomb counterterms. The latter is important to guarantee numerical 
stability and precision high enough to be able to give a quantitative statement 
concerning the (non)degeneracy of states. In fact, the convergence of the eigen¬ 
values turns out to be exponential in the number of integration points. It is 
plotted in Fig. 

The results are compared to equal time perturbation theory (ETPT). The 
spectrum is shown in Fig. ||, which is the most important of this article. A one- 
to-one comparison between the multiplets for a Bohr quantum number n = 2 
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Fig. 2. — Positronium spectrum in different Jz sectors. A^i = N 2 = 21. 


is depicted in Fig. ||. It is quite clear how to interpret these results. One reads 
off easily from the number of degenerate Jz levels the quantum number J of 
the complicated operator J^, consisting of the angular momentum operators, 
c/. the theory of the Poincare group in front form dynamics |^, Chapter 3.1]. 

It is worthwhile mentioning another non-trivial result of the formalism con¬ 
sidered in this article. The annihilation channel can be included straightfor¬ 
wardly, which is an intrinsic check of the effective theory. Surprisingly the 
instantaneous and the dynamic graphs act in different Jz sectors. An interest¬ 
ing feature concerning the annihilation channel is the calculation of the hyper- 
fine splitting, where this channel plays a major role and yields a well-known 
contribution. The coefficient of the hyperfine splitting is defined as 


Chf 


Mt-Ms 


1 

2 


2 

3 



The results for the two formalisms are 



a 

TT 





; this work 
;ETPT[0(a'‘)], 
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which are in reasonable agreement. 

Let us have a closer look at the results. One question listed above was, how 
good is the numerical evidence for degeneracies. One can plot the deviation of 
corresponding eigenvalues for Jz=0 and Jz=i multiplets with growing number 
of integration points, as was done in Fig. 12 of Ref. [||. A fit of the difference 
function to 

AM'^{N) = a — bexTp{{N — No)/c} . 

yields for the ground state triplet 1 ^Si: a = —(5.47± 0.95) x 10“®, b = (1.88± 
0.03) X 10“^, c = 4.03 ± 0.11. It seems thus quite justified to consider these 
states degenerate. 

The cutoff dependence of the eigenvalues was investigated, too. For the 
actual graphs, see Ref. . The role of the annihilation channel seems to be to 
stabilize the eigenvalues. For the triplet ground state we can fit the dependence 
on log K with 


f 3.90976 — 0.01858 log A + 0.00789 log^ A ;noannih. 

\ 3.91392- 0.00029 log A+ 0.00015 log^ A ;incl. annih. 


The decrease of the triplet with log A is suppressed by the inclusion of the 
annihilation channel by a factor of 60! Care is, however, to be taken, since 
the annihilation channel is absent in all QED bound states not built out of a 
particle and its antiparticle |Q. 



Fig. 3. — Comparison of multiplets for n=2: (a) results of the present work with 
a = 0.3, A = 1.0m/, = N 2 = 21; (b) equal-time perturbation theory (ETPT) up 

to order 0{a'^). 
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Fig. 4. — The positronium spectrum including the annihilation channel in the Jz=0 
sector. Parameters of the calculation: a = 0.3, A = 1.0 m/. The mass squared 
eigenvalues in units of the electron mass m? are shown as functions of the number 
of integration points N = N\ = N 2 - The triplet states, especially 1 are lifted up, 
the singlet mass eigenvalues are the same as without the annihilation channel. 


4. CONCLUSIONS 

To conclude, we review the results of the positronium theory in front form 
dynamics as presented in this article. 

The correct positronium spectrum is obtained in all Jz sectors. The Jz states 
form (degenerate) multiplets, which means that rotational invariance is restored 
in the solution of a Hamiltonian held theory in front form dynamics. Conse¬ 
quently, there is no need to diagonalize the complicated rotation operator 
because its eigenvalues can be read off by counting the number of degenerate 
eigenstates within the Jz multiplets. The annihilation channel can be intro¬ 
duced into the theory straightforwardly. Furthermore, the annihilation channel 
seems important for the cutoff behavior. From these facts, we can deduce that 
the applied underlying effective theory is correct. The computer code which 
was generated to solve for the eigenstates is applicable also to other effective 
Hamiltonians. The application to Wegner’s formalism ||] is work in progress. 
As an outlook one could think of plugging a running coupling constant into 
the code and to calculate the meson spectrum. Surely, this implies a deeper 
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Fig. 5. — The triplet ground state wavefunction for Jz= + 1 as a function of the 
longitudinal momentum fraction x and the transverse momentum k±, omitting the 
dependence on the angle (p. The calculation was done with a = 0.3, A = 1.0 irtf, Ni = 
41,A^2 = 11. Shown are: (a) (||)-component, (b) (|J,)-component. 


understanding of the numerical behavior of the occuring severe singularities. 
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